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Abstract 

Damping and thermal fluctuations have been introduced to collective normal 
modes of a magnetic system in recent modeling of dynamic thermal magneti- 
zation processes. The connection between this collective stochastic dynamics 
and physical local relaxation processes is investigated here. A system of two 
coupled magnetic grains embedded in two separate oscillating thermal baths is 
analyzed with no a priori assumptions except that of a Markovian process. It 
is shown explicitly that by eliminating the oscillating thermal bath variables, 
collective stochastic dynamics occurs in the normal modes of the magnetic 
system. The grain interactions cause local relaxation to be felt by the collec- 
tive system and the dynamic damping to reflect the system symmetry. This 
form of stochastic dynamics is in contrast to a common phenomenological ap- 
proach where a thermal field is added independently to the dynamic equations 
of each discretized cell or interacting grain. The dependence of this collective 
stochastic dynamics on the coupling strength of the magnetic grains and the 
relative local damping is discussed. 
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I. INTRODUCTION 



A new model has recently been developed [1], [2], [3] to study thermal noise and dynamic 
thermal reversal in interacting magnetic systems. In this approach, damping and thermal 
fluctuations are introduced to the independent normal modes of the magnetic system, cor- 
responding to the analogy of temperature defined by independent particles in an ideal gas. 
The damping term in the dynamic equations differ from that of Landau-Lifshitz [4] and, 
for even a single domain particle or film, reflects the asymmetry of the magnetic energy [3] . 
Generalization of the LLG equations to reflect the magnetic symmetry has been discussed 
in general [5] with specific analysis for the conduction electron relaxation process [6]. Col- 
lective normal mode processes have also been examined through analysis of relaxation to 
the complete spin- wave spectrum in thin films [7] . Here we derive the stochastic differential 
equations (SDE) specifically for the case of local damping in a system of interacting grains. 

Historically, stochastic differential equations have been developed by simply adding a 
thermal fluctuation field to the Landau-Lifshitz-Gilbert equations [8]. This approach has 
been widely utilized to analyze the role of thermal fluctuations, for example, in non-uniformly 
magnetized materials, such as a thin film, by discretizing the film and solving the coupled 
LLG equations for each cell with a statistically independent random field added to each 
cell [9], [10]. One argument in favor of this individual particle approach is that for physi- 
cally localized relaxation processes damping and thermal fluctuations can be conveniently 
introduced to individual particles as an effective field. 

However, from a collective normal mode point of view, even physically localized relaxation 
processes should give collective stochastic dynamics in the normal modes. Application of 
this approach to analyze thermal noise in a thin film appears to give better agreement with 
experiment than the LLG-Brown approach [11]. Here, collective stochastic dynamics are 
explicitly derived through a system-reservoir interaction model. We consider two coupled 
magnetic grains embedded in two separate thermal baths, focusing on small amplitude 
oscillations near equilibrium. No a priori assumption is made concerning the form of the 
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dynamic damping. This configuration provides a simplified picture for localized relaxation 
processes. We expand the analysis in [3] for a single grain utilizing the method of [13] 
to add a generalized thermal bath to the magnetization dynamics. The technique is to 
explicitly eliminate the oscillating bath variables to obtain a closed stochastic equation for 
the magnetic system. Under a Markovian approximation it can be shown that the magnetic 
system obeys collective stochastic dynamics in the form of damped harmonic oscillators 
in the normal modes. Thus, the original conjecture is verified that damping and additive 
thermal fluctuations should be added to the normal modes of a magnetic system, even if the 
physical relaxation processes are local. 

Section II introduces our model configuration. Section III obtains the dynamic equa- 
tions for the magnetic system by explicitly eliminating the bath variables. In Section IV 
Markovian and rotational wave approximations are utilized to obtain the collective stochas- 
tic dynamics for the independent normal modes of a magnetic system. Section V discusses 
the dependence of this collective stochastic dynamics on the magnetic interaction and a 
comparison to the individual particle picture is given. 

II. TWO INTERACTING GRAINS EMBEDDED IN DIFFERENT LOCAL 

THERMAL BATHS 

We consider two interacting nonidentical cubic magnetic grains of diameter D embedded 
in two different localized thermal bathes, as shown in Fig.l. Initially, neglecting the coupling 
to the thermal baths, the energy for the host magnetic grains, normalized by M^V, is: 

E m = M 2 s V\^h kl (l - m\ z ) + -h k2 (l - m\ z ) - h ex rtii • W? 2 (1) 
+(rft 1 ■ rrt 2 - 3m lx m 2x )} 

where rft 1 = A?x/M s and rft 2 = ~M 2 /M S are normalized magnetizations of the two grains. 
The normalized anisotropy fields are h kl = 2fT 1 /M s 2 and h k2 = 2K 2 /M 2 . The normalized 
exchange field, assuming an interaction strength A (acting uniformly between grain centers), 
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is h ex = 2A/D 2 M 2 (although here normalization is by M 2 rather than by K typical in 
micromagnetic simulations of recording media [14]). Normalization by M 2 V yields the 
magnetostatic term (using a dipole approximation) without any multiplicative constant. 
The magnetostatic and exchange interactions are assumed to be small compared to the 
anisotropy energies h k i,h k2 > 1 > h ex . With this assumption the equilibrium state is 
(Fig.l): 

77? i = e z , rrt 2 = e z (2) 

For small excitations around equilibrium (e 2 ), we only need consider second order variations 
so that: 

m iz = 1 - \( m \x + m iy) (3) 
m lz m 2z = 1 - ^{m\ x + m\ y + m\ x + m\ y ) 

The magnetic energy of the host grains, to second order, is: 

-(-/ifci + 1 - h ex )(m 2 lx + m\ y ) - - 
+ (1 - h ex )m ly m 2 y - (2 + h ex )m lx m 2x 



E m = M 2 V[--(-h kl + 1 - h ex ){m\ x + m\ y ) - -{-h k2 + 1 - h ex )(m 2 2x + m\ y ) (4) 



For each grain (j), we transform the two magnetization components orthogonal to the equi- 
librium direction ( mj X1 m^ y ) into (linearized) rotating magnetization components (a*,dj) 
(e.g. [12]): 

_ m 3X + im jy _ m jx - im jy 

In these coordinates, again in the lowest order quadratic variation, the magnetic energy (4) 
can be written as: 

E m = M 2 V[-(-h kl + 1 - h ex )aia[ - (-h k2 + 1 - h ex )a 2 a* 2 (6) 
+ 1/2(1 — h ex )(a\a 2 + a\a 2 — a\a 2 — a\a 2 ) 
— 1/2(2 + h ex )(aia 2 + a\a 2 + a\a 2 + a*^)] 



In order to consider localized relaxation processes the host magnetic grains are bilinearly 
coupled to two separate oscillating thermal baths [13]. This is a simplified model for local 
relaxation processes of the interacting magnetization system. A physical example would be 
relaxation by localized high moment Rare Earth impurities (e.g. see [15]). The total energy 
including magnetic energy, thermal bath energy and interaction energy is: 

Et = E m + E b + Ei (7) 

E m = E m (aia\, 02, a 2 ) 

M S V ^ M S V ^ 

&b — / , ^Ife0ifc0 lfc H > Ld2k02k0 2 k 

' k ' k 

M S V ^ M s y ^ 

^7 = 2^ dikihk^ + 6 lfc ai) H 2^ 92k{02ka 2 + b 2k a 2 ) 

^ k ^ k 

where bjk,b* k are the oscillating thermal bath variables. gik,92k represent the coupling 
strength of the magnetic systems and the thermal baths. Note that the bath terms are in 
the form of independent harmonic oscillators. Using the transformation (5), the interaction 
term Ei can be viewed as simply a Zeeman coupling with a thermal field. Here following 
[13], we assume that thermal bath is in equilibrium and the magnetic coupling is only a 
small perturbation to thermal baths. 



III. MODEL DYNAMIC EQUATIONS 

The standard procedure is to obtain a closed dynamic equation in the magnetization vari- 
ables by explicitly eliminating the thermal bath variables [13]. The Hamiltonian equations 
for the magnetic system are: 



ddj . 7 dE 

~~dt = ~ l M S V da* 

da) 7 dE 
= i- 



(8) 



dt M S V da 
which can be written in matrix form as: 
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E 9ikK k 
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— E 92khk 
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\a* 2 ) 




[a* 2 j 




T,92kb* 2k , 



(9) 



where: 



G = 



1 M S 
(10) 



( (-h kl + l-h ex ) (l + 2h ex )/2 3/2 ^ 

~(-h kl + 1 - h ex ) -3/2 -(l + 2*h ex )/2 

(l + 2h ex )/2 3/2 (-h k2 + 1 - h ex ) 

V -3/2 -(1 + 2 * h ex )/2 ~(-h k2 + 1 - ^ J 

is the matrix for gyromagnetic precession. 

In order to obtain a closed equation for magnetic system variables (a±, a*, a 2 , a^), we need 
to formally represent (bi k , b* lk1 b 2k , 6^) in terms of (a±, al, a 2 , a 2 ) in (9). We expect the closed 
equation for (a±, a*, a 2 , a 2 ) to be a stochastic differential equation under a Markov approxi- 
mation. However, before starting to represent (bi k , b* lk1 b 2k , b^ k ) in terms of (ai, a* X) a 2 , a 2 ), we 
notice that (9) has mixed times because of nondiagonal terms in the matrix (10). Thus, due 
to intergranular interactions, the distinct magnetic system time scales are not well presented 
in the rotating magnetization components (a±, a*, a 2 , a^). In the Markov approximation for 
stochastic dynamics, distinct system time scales must be separated from thermal bath time 
scales. So we first need to obtain explicitly system time scales [16]. This is done by normal 
mode analysis (e.g. [12]). 

The nondiagonal matrix (10) can be diagonalized into the following form: 



6 



G = UAU- 1 
( 



(11) 



U 



\ 



Un Vu U12 V 12 

Vll U U V12 U12 

U21 V 2 1 U 22 V 22 

\ V 2 1 U 2 1 V 22 U 2 2 J 



A 



^-U! ^ 

ui 

-uo 2 

y u 2 / 



Thus, the system gyromagnetic motion alone without a thermal bath or equivalently a 
relaxation mechanism can be written as: 

/ 



d 
di 



c 2 



«2 



c 2 



(12) 



-iwi 

iwi 

-iw 2 

iu 2 J 

un v u u 12 v 12 

Vu U U V12 U12 
«21 l»21 «22 ^22 
\ f 21 «21 ^22 «22 y 

lwhere (ci, c*, C2, c|) are the system normal modes and distinct system times scales can be 
determined as 1/cuiand l/u 2 . c^iand cj 2 are the magnetic system resonant frequencies. 
Now we represent total Hamiltonian (4) using normal mode coordinate: 



C2 
\C*2/ 



E 



M s V r , „. M S V ^ , M S V x ^ , ,. 
[wiCiCi + uj2C2C 2 \ H > uj lk bi k b lk H > ^2kb 2 kb 2k 

1 k ~ 1 k 

^ ci(gikVuhk + 9ikUnb* lk + g 2k v 2 ib 2k + g 2k u 2 ib* 2k ) ^ 



(13) 



7 



cJ(^lfcWll6lfc + fflfcUllftlfc + ^2fc«21&2k + g2kV2ib* 2k ) 
C2(gikVl2blk + gikUl2b* lk + g2kV22hk + g2kU22b* 2k ) 

\ c* 2 (gi k ui 2 b lk + gi k v 12 b* lk + g 2 kU 2 2b2k + g2kV22b* 2k ) J 
It should be pointed out here that diagonalizing the matrix G in (11) together with the re- 
quirement that the normal mode energy in form of E = ^^[ojiCic\ + uj 2 C2C2\ (13) determines 
the transformation matrix U uniquely. The dynamic equations for the magnetic system can 
be obtained from the Hamiltonian: 



so that: 



dt M s Vdc* 1 ; 



dc 

_I = -iu lCl - i ^2 (gikUubik + 9ikVub* lk + g 2 kU2ib 2 k + 9ikV 2 ib* 2k ) (15) 
k 

^ = -iu 2 c 2 - i ^2 (gikUuhk + gikVi2b\ k + gikU^k + g 2 kV 22 bl k ) 
k 



The dynamic equations for bath variables are: 



db ik . 7 dE 

- l T7J77^T ( 16 ) 



that is: 



dt M S V db* k 



<U ,fc = -iu lk b lk - i(^ifeWiiCi + gikVucl + SfifcWi 2 c 2 + #1*^12^) (17) 



dt 
db 2k 



-iUJ 2 b 2 k ~ %2jfc«2lCl + #2^21^ + 92kU22C2 + 92^22^) 



dt 

Equation (17) can be formally solved as: 

bi k (t) = hki^e'^ - tg lk [ dt'(u n ci + v n c{ + u 12 c 2 + v^De'^-^ (18) 

Jo 

b 2k (t) = b 2k (0)e- i ^ t - ig 2k [ dt'(u 2 i Cl + v 21 c\ + « 22 c 2 + v^d^e^^ 

Jo 

In order to simplify the calculation, in the following we assume the two thermal baths 
are identical except for a weighting factor in coupling strength to the magnetic system: 

92k = Pgik, uik = u 2k (19) 
Substituting (18) into (15), we obtain: 
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^ = -iu lCl - % 9ik[uub lk (0) + f3u 21 b 2k (0)]e- iuJkt (20) 

k 

k 

~ ^ 9ik / d£e~ Wkif ~ 1f) {u\ 1 c 1 + UuVucl + u u u 12 c 2 + u u v 12 c 2 ) 

k J ° 

- ^2 P 2g i k / dt'e'^'^iu^C! + u 2l v 2 ic\ + u 2 iu 22 c 2 + u 21 v 22 c* 2 ) 
k J ° 

+ Yl 9 lk / ^e^^^U^Ci + V u u n cl + + Ull«l2^) 

k J ° 

+ Y] I3 2 g 2 lk I dt'e-^-^iv^d + v 21 u 21 c{ + v 21 v 22 c 2 + v 21 u 22 c* 2 ) 
k Jo 

and similarly for c 2 . 

IV. STOCHASTIC DYNAMIC EQUATIONS 

We assume that the thermal bath is in thermal equilibrium with a continuum density of 
states: 

< Mo) >=< £4(0) >= o (21) 

< b ik (0)b* k ,(0) >= 5 kk m k 

where % = 1,2, n k is the energy of the kth component of the oscillating thermal bath, 
proportional to temperature if the thermal bath is in equilibrium. The summation over 
discrete energy levels can be approximated by a continuous integration: 

E 3ik - / du k D(u k )g 2 (u k ) (22) 

k 

and where occupation number in (21) is directly proportional to the temperature: n k ~ ksT. 

We need to single out the magnetic system time scales from the coupled total system 
time scales. This is done by the following transformation [13]: 
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ci(*) = Cl {t)e^\ c 2 (t) = c 2 {t)e^ (23) 
Substituting (22) and (23) into (20), we obtain: 

§ = fi(t)-fdu> k D(u> k )g 2 (u> k )fidtf 
+ /3 2 t^ 1 )ci(« / )e- i ( w *- wi )(*-*') + (Miifii + /5 2 M2if2i)cl(t')e- !(wfe+QJl)(( ^' ) e &li 



li 



(«11«12 + / 3 2 M2lM22)c 2 (t , ) e_i(Wfc ^ l)( '^' )eiK ^ 2) ' + («llUl2 + ^«2lf22)c2(t')e^ (Wt+ " l)(t_t ' ) e !(Ul+U2)i 

-(v 2 u + /3X 2 1 )? 1 (t')e i( ^ +Wl){ *-*' ) - (unuii + (3\ 2l v 22 )c\{t l )e l ^-^ t ~ t '">e 2i ^ t (2 
-(unui2 + /3 2 ^i^)^(Oe i(wfc ^ )( *~* ) e <(wi - W2) * - (unui2 + /3ViM22)c^(t , )e i( ^ +W2)( *~*' ) e i(wi+CJ2) ' 

where 

/i(f) = -i J] ffifc[«n6ifc(0) + Pu*ib2k(0)]e~ i{uk ~ Ul)t (25) 



) 



-< £ 9ik[vnb* lk (0) + Pv 21 b* 2k (0)}e^^ 



k 

is the thermal fluctuation field for the first mode. It satisfies: 



< h{t) >= (26) 

< fi(t)fm >= DMg'MnMiul, + P 2 u 2 2l )5(t - if) 

A similar expression for the second mode can be obtained. 

In the Markovian approximation [13], [16], we neglect memory for the slow variable 
(ci(i') — > Ci(t)). For long times t » l/(a>i — oj k ) the upper integral limit may be put to 
infinity and using: 

f°° 1 

/ dxe luJX = 7v5(uj) +iP.V.(~) (27) 
Jo v 

only the following terms remain: 
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f = h{t)e^ - J du k D{u k )g\u k ) dt> 
( M 2 1 + /5 2 M 2 i )c 1 (t)e- i ^-^)(*-*') 
(«n«i2 + /3 2 M 2 iM2 2 )c 2 (t)e- J ^-^)(*-*')e^^-^)* 
-(unwii + (3 2 v 2l v 22 )c ¥ 1 (t)e i ^-^ t - t ^e 2i ^ t 



(28) 



Because e 2 ^ 1 * is a fast oscillating term for t » l/^i, the — (unWn + 
(3 2 V2iV2 2 )c* 1 (t)e i ^- UJl ^ t - tn >e 2iuJlt term can be included into the fluctuating term. This is the 
rotational wave approximation [16]. If the gyromagnetic system is not degenerate (cu 1 ^ u 2 ) 
and the gyromagnetic rotating frequency gap is sufficiently large t » — lu 2 ), e 1 ^ 1- ^ 2 ^ 

is also a fast oscillating term and the (u 11 u 12 +P 2 U2iU22)c 2 (t)e~ i ^ k ~ U!1 ^ t ~ t ^e i ^ 1 ~ U! ^ t term also 
enters into the fluctuating terms. As discussed in section V, for a wide range of parameters, 
even for identical grains, nondegenerate gyromagnetic rotation is guaranteed. 
The following stochastic differential equation for c\ (t) is obtained: 



Weak interactions between the system and the thermal reservoir is usually assumed for 
a Markovian approximation and this corresponds to 771 << 1. The thermal fluctuation 
term(29) is: 



dci 
~dt 



F^t) - D{cu 1 )g 2 (u; 1 )(u 2 11 + /3 2 t&)c!(t) 



(29) 



where the damping term is defined as: 



t )l = D( Wl ) fl 2 ( Wl )« + ^ 1 ) 



(30) 




(31) 



k 



+ [ Vll bl k (0) + Pv 21 b; k (0)]e^ + ^} 

r poo 

~ / du k D(u k )g 2 (u k ) / ^( Mll «i2 + /9 2 w 2 iM 22 )c 2 (t)e- i ^- 1 )('-'') e ^— ^ 





+ / du k D{u k )g 2 {u k ) / di , ( UllUll +/3Vu22)cJ(t)e <(h '*- h ' l)(t - t ' ) e : 
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Notice that the last two terms are from the fast oscillating terms and have the magnitude 
proportional to rji. The thermal fluctuation has zero mean and, to leading order (neglecting 
r) 2 ), the variance is: 



< F 1 (t)F*(t') >= DMg'MnMiul, + (3 2 u 2 2l )5{t - t') 



(33) 



A similar equation for c 2 (t) can be obtained based upon the same conditions: 



?£ = F 2 (t)-T h &(t) 



(34) 



r] 2 = D{u 2 )g 2 (u 2 )((3 2 u 2 22 + u\ 2 ) 



F 2 (t) = -iY, 9ik{[u 12 b lk (0) + Pu^me-^-^ 1 



k 



+[v 12 b* lk (0) + (3v 22 b* 2 MV {Uk+u)2)t ) 



+higher order fast oscillating terms 

To summarize, the thermal fluctuation in the two modes have to leading order the following 
properties: 



The last condition in (35) gives uncorrelated thermal fluctuations in the two normal 
modes. It should be pointed out that in principle the fast oscillating terms in the thermal 
noise could result in correlation between Fi(t) and F 2 (t). However, those terms are of order 
i] 2 and can be neglected for a consistent Markovian approximation with weak interactions 
between the magnetic system and thermal reservoir. 

Transforming (29) and (34) into C\{t) and c 2 (t) coordinates, we obtain the following 
stochastic equations for the collective normal modes of the interacting magnetic system: 



<F 1 (t) >=0,<F 2 (t) >=0 



(35) 



< F 1 (t)F*(t') >= DMg^nMiu^ + (3 2 u 2 2l )5(t - t') 



< F 2 {t)F* 2 {t') >= D{u 2 )g 2 {u 2 )n{u 2 ){f3 2 u 2 22 + u\ 2 )5{t - t') 



< Fi{t)F 2 {t') >= 
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(36) 



-TT + V2C2 = ~iuJ 2 C2 + f2(t) 

at 

7]i = D(CU 1 )( ? 2 (CU 1 )«+/5 2 ^ 1 ) 

r] 2 = Diu 2 )g 2 iuj 2 )i(3 2 u 2 22 + u 2 l2 ) 
<fi(t) >=0,</ 2 (t) >=0 

< /i(*)/r(0 >= I>(wi)<T , (a;i)n(wi)(«? 1 + /Pt&M* - f ) 

< /i(0/ 2 *(O >= 



Thus, this analysis without any a priori assumptions has given stochastic differential 
equations for the two independent normal modes in the form of independent damped har- 
monic oscillators driven by uncorrelated thermal fluctuations. Note that the damping terms 
7]i iui) and r]2 iio 2 ) are directly proportional to thermal variance terms < fi(t)f*(t') >= D c \ 
and < f 2 {t)f2{t') >= D c2 respectively. The ratio of fluctuation magnitudes in two modes 
are: 



For a magnetic system that satisfies a canonical distribution around equilibrium, the mag- 
nitude of the damping coefficients 771 and 772 are related to the thermal fluctuations terms 
D c i and D c2 through the fluctuation-dissipation condition [13]: 



D c2 



D c i 



Diui)g 2 iui)niui)iu\ l + p 2 u 2 21 ) 
D(u 2 )g 2 (uj 2 )n(uj 2 )(P 2 ul 2 + u\ 2 ) 



(37) 



Vi = 



DciM s V 
lk B T 

D c2 M s V 
lk B T 



U 2 = OiQ 2 U 2 



■ui = a Q iUi 



(38) 



where rji and 7] 2 are damping rate. 



V. DISCUSSION 



We have shown here (36) that local physical relaxation mechanisms give collective 
stochastic dynamics in the system normal modes for interacting magnetic grains or con- 
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tinuum discretization cells. For a given system the specific damping terms can be evaluated 
using (38). These results can be generalized for any system of N interacting magnetic sub- 
units. The general operative approach is to first diagonalize the (2N x 2N) matrix of the 
gyromagnetic precession (near equilibrium) without any damping terms to find the collective 
modes and their eigenvalues or resonance frequencies. Following that, using expanded forms 
of (38), damping and thermal fields are added to each of the collective mode dynamics. 
The total energy in each mode is ksT and the fluctuation-dissipation theorem relates the 
variance of the thermal term to the damping. As will be shown in the specific examples be- 
low, there is no contradiction between local relaxation mechanisms and collective stochastic 
dynamics. The interactions between the elements allow damping in one cell to be felt by all 
the others. The interactions also cause the dynamics to reflect the overall sample geometry 
or anisotropy. 

Here we evaluate some examples for this two-grain system. The most general case of the 
dependence of the collective stochastic dynamics on interaction strength can be obtained 
from (11). We begin with the results for a system of identical grains with no exchange 
(h ex = 0) . The mode resonance frequencies are: 



In Fig. 2a these normalized frequencies are potted versus interaction strength M S /H K1 . 
With weak or vanishing magnetostatic interactions the frequencies are simply that of inde- 
pendent grains with cu 1 = ^H K1 and uj 2 = ~{H K 2 = ~{H kl . With increasing interaction both 
resonance frequencies decrease. The modes in zero field are asymmetric coherent rotation 
(mode 1) and asymmetric fanning (mode 2) [17]. The lowest frequency mode is almost coher- 
ent, because there is less magnetostatic energy to rotate into the particle axis direction (x: 
Fig.l). As the magnetostatic interaction is increased the energy barrier (or quadratic curva- 
ture) decreases and at M s /H K i ps 0.33 the frequency vanishes and the equilibrium direction 




(39) 
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moves from the "z" direction to be along the line joining the particles (the "x" direction). 
For the higher energy mode the frequency uo 2 also decreases with interaction strength and 
eventually the energy barrier will vanish (002 — > 0), but at a higher M s /Hki = 0.5 due to 
the larger magnetostatic energy of that mode. 

In Fig. 2b normalized resonance frequencies are plotted versus exchange h ex for the case 
M s /Hki = 0.09 and Hk2 = 1.95Hki- At h ex = the resonance frequencies are almost the 
ratio of the Hk values. The quasi-coherent mode frequency (coi) hardly varies with exchange, 
as expected, since coherent motion does not involve exchange energy. The incoherent mode 
frequency 002 increases rapidly with exchange because of the increased exchange energy of 
that non-uniform mode. We emphasize that the stochastic dynamic modeling in this paper 
is for 0J2 7^ cc>i(and sufficiently different). However, the plots shown here indicate that 
for almost all cases of interest, this condition holds, even for identical particles with finite 
coupling. 

We now explore the variation of the mode damping parameters with magnetostatic in- 
teraction for h ex = 0, as in Fig.2a and with slightly differing anisotropy (H K2 = l.OOlif^i). 
As can be shown using (11), for M s /H K i = , we have ttn = u 22 = 1, ^11 = v 22 = 0, 
U21 = v 2 i = 0, U12 = V\2 = 0. The normal modes are just those of the individual particles: 
(ci = ai, c 2 = a 2 ) and the damping terms (36) are just the local damping for the individual 
grains: 

77! = L>(^i)<?V) (40) 
772 = f3 2 D(u 2 )g 2 M 

In Fig. 3a we plot normalized damping terms D(ui)g 2 (ui)n(LJi), 112/ D(u 2 )g 2 (002)11,(002) 
versus magnetostatic interaction M s /Hxi for the case of /3 — 1. For M s /Hx\ = 0, the 
damping terms are just those of the independent grains. As seen in (38) the damping terms 
are proportional to the mode frequencies and that variation is seen in Fig. 3a. It is noteworthy 
that for identical damping mechanisms and a finite grain interaction, the collective modes 
are not identical. 
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The case for (3 = 0, corresponding to relaxation only in one cell, is shown in Fig.3b. 
In this example with no interactions (M S /H K1 = and h ex = 0), i] 2 = corresponds to 
no thermal bath coupling for grain 2. For a very slight non-uniformity {Hki = 1-OOlf/xi), 
as M s /Hki increases, the normal modes become a combination of the individual particles 
and for (H K2 — H K1 ) /H K1 = 0.001, this occurs when the magnetostatic field surpasses the 
difference of the anisotropy fields (M s /Hki ~ 0.001) . For M s /Hki > 0.001 the normalized 
damping terms become almost equal as the coupling dominates (r]i ~ r\ 2 ~ 0.5) and then 
decrease following the frequency dependence in (38). 

In Fig. 4 the damping terms are plotted versus exchange for a fixed magnetostatic interac- 
tion M S /H K1 = 0.09 and H K2 = 1.95H Klt as in Fig. 2b. In this case f3 2 = 0.25, corrsponding 
to a weaker relaxation in grain 2 compared to grain 1. The initial values for h ex = re- 
flect the frequency dependence of the damping, between the values as seen in Figs. 3a, 3b. 
With increasing exchange the damping of mode 1 hardly varies and the damping of mode 2 
increases rapidly, following the results of Fig. 2b and the accompanying discussion. 

As a comparison of the results of this analysis with the conventional LLG-Brown ap- 
proach, following [3], we compare the spectra of thermal white noise excitation of this sys- 
tem. We take the case of damping only in mode 1 (@ — 0). The analysis is in Appendix I. For 
M s /Hki = 0.25, H K2 = Hki, and h ex = 0.5, following Fig. 4 for the collective analysis, we 
use r] 2 /r]i = 1.5. For the LLG-Brown equations we use ct\ = 0.1, a 2 = 0. Both curves have 
been normalized by matching the first resonance peak and linewidth (by adjusting r\ 2 ja\). 
Both curves exhibit the same resonance frequencies for the two modes; the resonance fre- 
quencies are (to first order) independent of damping. The spectral shapes differ, however. 
As in [3], the low frequency PSD of the LLG approach is 3-4 dB above that of the collective 
result (greater than 3 dB because of the proximity of the two resonance frequencies). The 
second mode peak is lower and broader for the collective model. 

In this work the problem of introducing damping and thermal fluctuations for an inter- 
acting magnetic system is addressed using a physical model of system-reservoir interactions. 
No a priori assumption is made concerning the form of the dynamic damping. This ap- 
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proach is quite different from the LLG-Brown approach [8], where a dynamic damping is 
assumed in the LLG format for each individual grain or discretization cell. In the LLG- 
Brown approach, the underlying physical processes of damping and thermal fluctuations 
are not explicitly considered. There have been two papers [18], [19] that have attempted 
to justify the LLG-Brown model for interacting magnetic units. In [18] it is argued that 
the LLG-Brown approach gives dynamic equations that can be cast in the form of general- 
ized Newtonian dynamics (specifically in the form of an RLC circuit of coupled oscillators). 
However, this argument is only inferential and is not derived using any basic physics model. 
The result is equivalent to assuming the application of LLG-Brown to individual grains or 
discretized cells of a continuous medium, as is also the essence of [19]. These analyses do not 
derive a stochastic differential equation with damping and thermal fluctuating terms from 
system-reservoir interactions. 

The arguments in [18], [19] for independent thermal fluctuations for independent par- 
ticles in an interacting magnetic system is, in fact, a thermodynamic consistent condition 
required for localized LLG damping for individual magnetic units. However, in this paper, 
we have explicitly constructed a system-reservoir model without any a priori assumption of 
the damping format or thermal fluctuations. The explicit elimination of thermal reservoir 
variables in the model gives stochastic dynamic equations of a damped harmonic oscillator 
driven by thermal noise only in the collective normal modes. As we show in Fig. 5, different 
thermal noise power spectra are obtained. The analysis presented here verifies that damping 
and additive fluctuations must be added to the normal modes of an interacting magnetic 
system, even if the physical relaxation processes are local. 

VI. CONCLUSION 

A fundamental analysis of the stochastic dynamic equations for local coupling to a ther- 
mal bath has been performed for a system of two interacting grains. Here no a priori 
assumptions about the form of the dynamic damping term have been assumed. The results 
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are in the form of damped harmonic oscillators driven by thermal fields in the collective 
modes of the system, a result that previously was derived rigorously for a single grain with 
anisotropic energy variations about equilibrium. The LLG-Brown formalism in which a 
thermal field is added to each grain or discretization cell of a continuous medium is shown 
not to apply. Noise power spectra are evaluated and it is shown that the two resonance 
frequencies are broadened by the collective damping terms. Localized relaxation is felt by 
the system collective stochastic dynamics due to the intergranular interactions. The form of 
the dynamic damping term, also due to the interactions, reflects the overall system symme- 
try. For approximately identical grains these damping constants are about equal even if the 
physical damping occurs in only one of the grains. Noise power spectra are shown which give 
significant differences for the two models, thus providing guidance to future experimental 
analysis. 
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VII. APPENDIX I 

The stochastic differential equations (36) for two independent normal modes can be 
written as: 

^ + »7ici = -iwici + /i(t) (41) 
dc2 

-TT + V2C2 = ~iuJ 2 C 2 + f 2 (t) 

at 

< /i(*)/r(0 >= ma 

< / 2 (t)/ 2 *(0 >= 2D c2 

where the ratio of fluctuation magnitudes is: 

D cl D{u l )g 2 {u l )n{uj l ){ul l + f3 2 u 2 21 ) 



D c2 D(uj 2 )g 2 (uj 2 )n{uj 2 )((3 2 ul 2 + u\ 2 ) 



(42) 
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For magnetic systems that satisfy a canonical distribution around equilibrium, the mag- 
nitude of damping coefficients r\\ and r\ 2 are related to the thermal fluctuations D cl and D c2 
through the fluctuation-dissipation condition: 



D cl M s V 



lk B T 

D c2 M s V 
r] 2 = — , „ uj 2 = a Q2 u 2 



(43) 



The spectra of the normal modes can be calculated based upon (41). Here the calculation 
is done in a non-dimensional format and the frequency is normalized by 7M5: uj = u/^M s . 
The spectral densities for two modes are: 

1 2 



Si(cu) = < c\(uj)ci(uj) >- 
S 2 (uj) = < c* 2 (uj)c 2 (uj) > 



(-JU+JW1 - OLQiUJi) 

1 



(-juj+ju 2 - a Q1 u 2 ) 



2a Q1 k B T 

M*V 
2a Q2 k B T 



The magnetization can be represented by the normal modes as: 



^ mix ^ 



m 2x 
\m 2y J 



PU 



( c \ 



\c 2 ; 



where: 



(44) 



(45) 



P = 



( 1/V2 1/V2 ^ 

l/y/2 l/y/2 

-i/y/2 i/y/2 
\ -i/y/2 i/y/2 ) 

The spectral density function for the magnetization is defined as: 



(46) 



S{uj) =< rft(uj) + rft(uj) > 



(47) 
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where + denotes conjugate transpose and rft = [mi x ,m 2x mi y ,m 2y \. The magnetization 
spectral density is related to the spectral density of the normal modes through the transfor- 
mation matrix: 



TR = U + P + PU 

Thus, the magnetization spectral density is: 

S(u) = [TR 11 S 1 (u) + TR 22 S 1 (u) + TR 33 S 2 (u) + TR u S 2 (u)\ 
a Q1 (TR u + TR 22 ) a Q2 (TR 33 + TR U 



2k B T 
M 2 S V 



(48) 



(49) 



K-jUJ+jUJ! + a Q1 Ldi)\ +juj 2 + «Q2^l) 

For the LLG model, damping and thermal fluctuations are added to individual particles. 
For each individual particle: 



di 



= — 7A? x ~H — 0:7 



A$ x (Jvt x it) 



(50) 



dt '" M s 

where it includes the effective field —dE/dJt and the thermal fluctuating fields lt T . lt T ' 
satisfies the fluctuation-dissipation condition: < Hj(t) >= and < Hj(t)Hj(t') >= 
{2ak,BT /^M S V) 5(t — t')5ij. The normalized LLG equations to leading order are: 



d_ 

dt 







^ mix ^ 




/ 




m 2x 


= 1 M S A 


m 2x 






H 2y 






+ 7 




m ly 




m ly 






~ H lx 


\m 2y j 




\m 2y j 




V 


-K j 



(51) 



where: 



A 



ai(-h k i + 1 - h ex ) «i(2 + /i ex ) -h kl + 1 - h ex 

a 2 (2 + h ex ) a 2 (-h k2 + 1 - h ex ) -(1 - h ex ) 
-(-hki + 1 - h^) -( 2 + h ex ) oti(-h kl + 1 - h ex ) 



-(1 - h ex ) 
-h k2 + 1 — h ex 
-a^l - h ex ) 



\ -(2 + h e 



(52) 



~(-h k 2 + 1 - h ex ) -a 2 (l - h ex ) oi 2 (-h k2 + 1 - h ex ) J 
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and ot\ and a 2 are different damping parameters for the two particles. The nonzero correla- 
tions between fluctuation terms are: 



<<(t)HUt')>= 2 -^S(t-t') (53) 



<<(;)<(;') >=^_^-o 

< H T lx (t)Hl(t') >= 2 -^r5(t-t>) 



2a 2 k B T 



lM s V 
2a x k B T 

lM s V 

Using (51) , (52), we can calculate the nondimesionalized correlation matrix as: 



< H T 2x (t)H T 2x (t>) >= 2 -^r5(t-t>) 



C ^) = ,- T l yi I \ ' ^4 IT (54) 

(juI-A+) (-jujI-A) 

where / is the unit matrix. The spectral density can be obtained from the correlation matrix 
and the fluctuation magnitudes (53): 



2k B T 

S(uj) = [aiCoru + a 2 Cor 22 + ctiCor^ + a 2 Cor u ] (55) 
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